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Abstract. Galaxy Clusters (GCs) are the largest reservoirs of both dark matter and cosmic 
rays (CRs). Dark matter self-annihilation can lead to a high luminosity in gamma rays and 
neutrinos, enhanced by a strong degree of clustering in dark matter substructures. Hadronic 
CR interactions can also lead to a high luminosity in gamma rays and neutrinos, enhanced by 
the confinement of CRs from cluster accretion/merger shocks and active galactic nuclei. We 
show that IceCube/KM3Net observations of high-energy neutrinos can probe the nature of 
GCs and the separate dark matter and CR emission processes, taking into account how the 
results depend on the still-substantial uncertainties. Neutrino observations are relevant at 
high energies, especially at > 10 TeV. Our results should be useful for improving experimen- 
tal searches for high-energy neutrino emission. Neutrino telescopes are sensitive to extended 
sources formed by dark matter substructures and CRs distributed over large scales. Recent 
observations by Fermi and imaging atmospheric Cherenkov telescopes have placed interest- 
ing constraints on the gamma-ray emission from GCs. We also provide calculations of the 
gamma-ray fluxes, taking into account electromagnetic cascades inside GCs, which can be 
important for injections at sufficiently high energies. This also allows us to extend previous 
gamma-ray constraints to very high dark matter masses and significant CR injections at very 
high energies. Using both neutrinos and gamma rays, which can lead to comparable con- 
straints, will allow more complete understandings of GCs. Neutrinos are essential for dark 
matter annihilation channels like XX ~~ ?" / i+ / Lt_ ) where the neutrino signals are larger than the 
gamma-ray signals, and for hadronic instead of electronic CRs, because only the first leads 
to neutrinos. Our results suggest that the multi-messenger observations of GCs will be able 
to give useful constraints on specific models of dark matter and CRs. 



Contents 



1 Introduction 



1 



2 Galaxy clusters as high-energy particle sources 

2.1 Annihilating dark matter 

2.2 Cosmic rays 

2.3 Neutrino and gamma-ray emissions 



3 

3 
6 

10 



3 Neutrino and gamma-ray constraints 

3.1 Constraints on annihilating dark matter 

3.2 Constraints on cosmic rays 



12 

13 
16 



4 Summary and discussion 



19 



1 Introduction 

Galaxy clusters (GCs) are the largest virialized objects in the universe [see reviews 1, 2]. 
Their mass (~ 10 14 — 10 15 M@) consists mainly of dark matter. Many galaxies (typically 
~ 50 — 1000), including a couple of active galactic nuclei (AGN), are contained in GCs [e.g., 
3]. GCs are also sources of X-rays that are mainly produced via thermal bremsstrahlung 
emission from the hot intracluster medium (ICM). 

The existence of dark matter is known by its gravitational effects [4], but its particle 
properties are unknown. Many dark matter models are motivated by elementary particle 
physics, and the most popular candidate is weakly interacting massive particles (WIMPs) [5]. 
In supersymmetric models [6], the lightest supersymmetric particle can be stable, which is 
achieved by a discrete symmetry such as R-parity, although meta-stable dark matter is also 
allowed as long as the lifetime is sufficiently longer than the age of the universe. WIMP 
annihilation cross sections are consistent with the values required for the correct thermal 
relic density. An important way to look for dark matter signatures is by observing indirect 
gamma-ray and neutrinos signals from self-annihilation. These may be seen from various 
astrophysical objects, including the Galactic center, nearby galaxies, dwarf galaxies, and 
isotropic diffuse backgrounds [e.g., 7, and references therein]. The annihilation signal is 
proportional to the density squared, and not only the smooth dark matter halo but also 
substructures are relevant [8-10]. 

GCs can also be attractive sources to look for annihilation signatures of dark matter [e.g., 
11—13]. GCs are the largest gravitationally-bound objects in the universe, and are still 
forming today. Recent N-body simulations have indicated that small, bound structures in 
larger halos are less affected by tidal stripping than in smaller halos, and a large population 
of dark matter substructures is present [9, 10]. This may be the case especially at the cluster 
outskirts [14]. If the smallest halo mass is order of the Earth mass, as suggested by arguments 
of the kinetic decoupling of WIMPs in the early universe [e.g., 15], then the annihilation signal 
can be boosted by a factor of ~ 1000 for massive GCs [14, 16-21]. This strongly motivates 
searches for annihilation signals from GCs. Recently, Fermi has started to give stringent 
constraints on the annihilation cross section for low-mass WIMPs [19-25]. As indicated in 
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Table 1. Comparison among dwarf galaxies, normal galaxies and GCs. For dark matter annihilation, 
GCs have an advantage in the boost factor since many substructures may survive. The boost factor is 
estimated based on Ref. [14]. For CRs, GCs have an advantage that higher-energy CRs can be trapped 
because of their larger size. The magnetic field is assumed to be 1 fiG for simplicity. Although the 
numbers are crude, this table demonstrates the importance of GCs for probing high-energy emissions. 



Source 


mass [Mq] 


Boost factor 


size [pc] 


E = eBr [eV] 


Dwarf Galaxies 


~ 10 Y 


~ 1 


~ 500 


- 5 x 10 17 


Normal Galaxies 




- 80 


- 2 x 10 4 


- 2 x 10 19 


Clusters of Galaxies 


~10 15 


- 1000 


- 10 e 


~ 10 21 



Table 1, one of the key points in gamma-ray and neutrino searches from GCs is the large 
boost factor allowed by the existence of substructures. 

GCs are also the largest cosmic-ray (CR) reservoirs in the universe. Various sources 
have been considered in the literature, including AGN [26-28], supernova remnants [29] 
and ram-pressure striping of infalling galaxies. Among many possibilities, AGN jets and 
shocks associated with the large scale structure formation seem the most promising [30- 
32, 34]. In hierarchical scenarios of cosmological structure formation, GCs should continue 
to grow through merging and accretion of dark matter and baryonic gas, thereby generating 
powerful shock waves on ~ 1 Mpc scales [33, 35-42]. At CG outskirts, around the virial 
radius, filaments and voids, shocks with relatively high Mach numbers are expected due to 
accretion shocks and minor merger events, and efficient acceleration of high-energy particles 
is possible [35, 37, 41, 42]. Low Mach number shocks, involving merger and flow shocks, 
occur inside the virial radius, which is important for thermalization [32, 39, 41, 42]. 

These CRs interact with the ICM via hadronuclear processes such as the pp reaction, 
leading to gamma rays and neutrinos [26, 31]. CRs may be accelerated above ~ 10 17 eV [28] 
and even up to ultra-high energies (UHE; > 10 18 ' 5 eV) [43] , and they might contribute to the 
observed CR flux at the Earth. CRs with ~ 10 18 — 10 19 eV interact with cosmic microwave 
background (CMB) photons via the Bethe-Heitler process, which leads to GeV-TeV gamma 
rays via the electromagnetic cascade in GCs [44, 45]. Photomeson production, which leads 
to gamma-ray and neutrino production, may also happen if even higher-energy CRs exist 
in GCs [45-47]. Predictions for the associated nonthermal gamma-ray radiation, including 
leptonic emissions, have been discussed by a number of authors. As indicated in Table 1, 
one of the key points is that higher-energy CRs can be trapped in GCs because GCs have 
~ 0.1 — 1 /iG magnetic fields over ~ 1 Mpc . 

Recent gamma-ray observations have started to allow significant probes of dark mat- 
ter [19, 21, 22, 48-50] and CRs [18, 48, 50-52] from GeV to multi-TeV energies, as expected 
from papers from the pie-Fermi era [53-55]. In the GeV range, Fermi has given interesting 
limits on models of dark matter and CRs. In the TeV range, imaging atmospheric Cherenkov 
telescopes (IACTs) such as H.E.S.S., MAGIC, and VERITAS play an active part in detecting 
very- high-energy (VHE) gamma-ray sources [56]. 

Neutrino observations should also be relevant at high energies, especially at > 10 TeV [57- 
59]. Neutrinos are special since they provide complementary probes of heavy dark matter 

1 Note that the proton maximum energy achieved by accretion shocks is much smaller than the value in 
Table 1 since the shock velocity is much smaller than the light velocity and the magnetic field would be weaker 
at larger scales [e.g., 43, 44]. 
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and clues to confinement properties of high-energy CRs. The IceCube detector at the South 
Pole was completed, and the comparable KM3Net detector in the Mediterranean sea is be- 
ing planned. Neutrinos have become one of the important and powerful messengers, along 
with gamma rays [e.g., 60], but there are only a few works about neutrinos from GCs in 
the contexts of both dark matter and CRs [23, 26, 28, 31] despite many on the gamma-ray 
emissions. Both gamma rays and neutrinos are essential for understanding high-energy pro- 
cesses in GCs. This paper focuses on the power of current and future high-energy neutrino 
observations. 

Confirmations and tests of predicted signals are possible by multi-messenger detectors. 
In particular, neutrinos provide a complementary probe even if gamma rays are more easily 
detected. Neutrinos can probe hadronic signals. For CRs, gamma rays may come from 
primary electrons, but neutrinos require protons or nuclei. For dark matter, neutrinos are 
produced by channels involving decay into heavy leptons, quarks, massive bosons, Higgs, 
and neutrinos themselves. Neutrino detectors have better sensitivities especially at higher 
energies. At such high energies, gamma rays are attenuated or cascaded whereas neutrinos 
can leave a source without attenuation, which allows us to reveal injection energies directly. 
Neutrino detectors typically have degree-scale angular resolution, and this is comparable to 
the size of nearby GCs. While the large angular size means that the IACT sensitivity is 
worse than its point source sensitivity, neutrino detectors have some advantages in observing 
such extended sources. 

The organization of this paper is as follows. In Section 2, we overview gamma-ray 
and neutrino production processes in GCs, and describe basic concepts. In Section 3, the 
results are shown. Summary and discussion are found in Section 4. We adopt Hq = lOOh = 
70.2 km s' 1 Mpc~\ fl dm = 0.229, tt m = 0.275, and A = 0.725. 

2 Galaxy clusters as high-energy particle sources 
2.1 Annihilating dark matter 

Many dark matter models have principal annihilation channels leading to quarks, W/Z gauge 
bosons, and Higgs. One may consider leptonic or semi-leptonic decays. For hadronic decays, 
many mesons are produced due to hadronization, leading to neutrinos, gamma rays and 
electrons/positrons via decay of mesons, and the gamma-ray spectrum from final states is 
near-universal though details depend on particle physics models [e.g., 5, 18, 21]. In addi- 
tion, one may have internal bremsstrahlung, i.e., direct emission from a virtual, charged, 
exchanged particle [e.g., 18, 61, 62, and references therein]. The internal bremsstrahlung is 
particularly relevant if lowest-order annihilation processes are suppressed e.g., via a helicity 
or loop suppression. It also allows electromagnetic emissions even for the vv channel in the 
tree level [61, 62]. 

In this work, for demonstration purposes, we use a representative annihilation channel, 
XX ~~ bb, where we run PYTHIA to evaluate the spectrum, dS/dE [63]. Examples of elec- 
tromagnetic and neutrino spectra are shown in Figures 1 and 2. We also consider another 
representative annihilation channel, XX A* A* - > motivated by leptophilic models. In those 
models, dark matter exclusively annihilates into a lepton and anti-lepton, which was exten- 
sively studied in the context of the PAMELA and Fermi excess [see reviews 7, and references 
therein]. Such channels may be realized in some models such as a Sommerfeld-enhanced 
model including light force-scale carrier [64]. For demonstrative purposes, we simply use the 
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Figure 1. Injection spectra of gamma 

rays (thick curves) and electron/positrons (thin 
curves) from dark matter annihilation (m x c 2 = 
5 TeV) for xx -> XX W + W~ , XX 

and XX ~ ^ 
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Figure 2. Injection spectra of neutrinos (sum of 
all flavors) from dark matter annihilation (m x c 2 = 
5 TeV) for xx -> P + ^ , XX -> XX -> ** 

and XX ~* bb. 



well-known formula of three-body decay of muons. The spectra are harder than those of 
channels involving hadronization. 

For comparison, spectra from other channels {\X ~^ W + W~ and XX ~> tt) are a l so 
shown. For the W + W~ channel, the weak boson decays into a neutrino and lepton and can 
also be hadronic, while the top quark mainly decays into a weak boson and bottom quark for 
the tt channel. In both the cases, the spectra are harder than that of the XX ~^ bb channel 
but softer than that of the fJ- + ^~ channel, so the resulting constraints would lie in between 
the bb case and the /J- + fJ>~ case. Also, we do not consider the internal bremsstrahlung since 
its relevance is model-dependent. 

Another particle physics factor is the self-annihilation cross section, < av > xx , which 
we attempt to constrain in this work. If dark matter is a thermal relic, the required cross 
section is ~ (2 — 5) x 10 -26 cm 3 s , depending on mass [65]. However, larger or smaller cross 
sections are allowed if dark matter is nonthermally produced in the early universe. The cross 
section is unknown, but neutrino and gamma-ray limits give useful constraints on < av > X x- 

In order to calculate indirect signals of annihilation dark matter, we need to evaluate 
the astrophysical factor as well as the particle physics factor. The astrophysical factor is 
determined by the dark matter distribution, which is not certain at present. The smooth 
dark matter profile is often described as [1, 5, 66] 

Pv(f) = Ps 7 71 1 (2.1) 

X (r/r s ) Q2 [l + (r/r s ) ai ] (Q3 " a2)/ai 

where p s and r s = r^ c jc s are the scale density and scale radius, respectively. Here is the 
virial radius such that the enclosed dark matter mass is Ma c = (47r/3)(A cy o c )r 3 ^ c , and c s is the 
mass-concentration parameter for given cluster mass. Note that A c is the overdensity relative 
to the critical density p c [1], and we use A c = 200 throughout this paper. Based on recent 
cosmological simulations, we also adopt the following mass-concentration parameter [e.g., 
67], 

7.85 / M 200 



C - (M2O ^ )= (TT^l 2xlO^M ' ^ 
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In Eq. (2.1), a\ is the parameter that determines the shape at r s , a% is the inner slope, and 
03 is the outer slope. The Navarro- Frenk- White (NFW) profile is given by ol\ = 1, ol% = 1 
and 03 = 3 [68], but different values are taken for other profiles such as Kravtsov [69] and 
Moore [70] profiles. Recent simulations seem more consistent with shallower profiles such as 
the Einasto profile [71], which also seems to be supported by some observations. On the other 
hand, it is natural to expect that the dark matter profile is steepened by baryon dissipation 
so that annihilation signals are enhanced [21, 72]. For reference, we consider the NFW profile 
throughout this work, but one should keep in mind that there is uncertainty in the smooth 
profile. 

Once the particle physics factor and astrophysical factor are evaluated, the intensity of 
gamma rays and neutrinos coming from the smooth dark matter halo (in the unit of particles 
per energy, area, time, and solid angle) can be written as 

™ (2 - 3) 



where r = yl 2 + (c?a©) 2 , W = y r A c ~~ (c^O) 2 , and O is the angular size measured from 
the direction to the CG center. The normalization of the spectrum is given by £ j dE E(dS/dE) 
i x c 



2m x c 2 . By performing integration over O, the total flux is obtained as 



? sm 1 < aV >XX dS 



where we have introduced the flux multiplier for given cluster mass [73], 

VjdV&r) 

s " iz) = rj^wF (2 ' 5) 

Here, V is the volume of the CG. The uncertainty in the smooth profile is included in this 
flux multiplier. 

In addition to the smooth dark matter distribution, there are contributions from dark 
matter substructures. To see interesting signals of annihilating heavy dark matter from GCs, 
we typically need significant boosts due to substructures, so that the smooth halo component 
is overwhelmed. The expected boost is especially relevant for larger halos, including GCs. 
The substructure model presented in Ref. [14] gives 

^(e ) ^(M 2 oo)Fr jln 1 (17)e2 + ( ; M/4)2 (2.6) 

where ©200 is the angular size corresponding to r2oo- We use the following boost factor, 

/ M \ °- 39 

&(M 2 oo) = 1.6xlO- 3 (^J , (2.7) 

where the minimum halo size is assumed to be M m \ n = 1O~ 6 M0. The boost factor also 
depends on the minimum halo scale as b oc M~j^ 2 [14]. Note that the total flux from dark 
matter annihilation is written as Fe = (1 + b)F E m . 

The gamma-ray and neutrino angular profiles are relevant for the detectability. To 
demonstrate this, we define the following quantity, 

J(< 0) = 2vr / dQ' I E (&). (2.8) 
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Figure 3. The flux (within 0) divided by 0, J(< 0), as a function of 0, which represents the shape 
of the signal to the square root of the background at given energy. We assume M min = 1O _6 M , 
where the contribution from substructures (thick curve) is largely dominant over the contribution 
from the smooth halo (dot-dashed curve). For the dark matter profile, the parameter set for Virgo 
is used [17]. The typical angular resolution of IceCube (0 = 1 deg) and CTA (0 = 0.05 deg) are 
indicated by the dashed vertical line and the dotted vertical line, respectively. The detectability is 
optimized at a few degrees in the presence of substructures. 



The results are shown in Figure 3, where substructures are included. Note that J(< Q)/Fe 
represents the normalized flux of the signal. It is further divided by © since the background 
flux is proportional to J7 = 2-/r(l — cos©) ~ 7r@ 2 , and we consider the signal relative to 
background uncertainties, i.e., the square root of the near-isotropic background intensity. 
Figure 3 indicates how searches for the extended emission are crucial in the case of GCs. 
The contribution from the smooth halo is concentrated in the central region, while that from 
the subhalos is much more extended. For nearby clusters such as Virgo, in the presence of 
substructures, the apparent size is ~ 1 degree, which is comparable to the angular resolution 
of IceCube. The optimization of the angular scale improves the detectability, though the best 
angular size depends on details of substructure models. For example, in the case described 
in Eq. (2.6), the signal from Virgo is maximized at ~ 2 — 3 degrees. 

2.2 Cosmic rays 

For astrophysical emissions, we focus on gamma-ray and neutrino production from the pp 
reaction, i.e., interactions with CRs and target ICM gas. This is enough as long as the 
CR spectral index is steep enough or the CR maximum energy is sufficiently less than ~ 
10 19 eV [28, 45]. The target gas density is related to the electron density, which is given by 
X-ray observations. The j3 model gives [1] 




(2.9) 
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where n co is the core ICM density, r co is the core radius and /3 describes the slope. The 
target nucleon density is determined by tin = n e /(l — 0.5xHe)j where x^e = 0.24 is the 
Helium fraction [e.g., 31, 54]. 

The total CR energy, £ cr , is unknown, so it is introduced as a free parameter that 
we attempt to constrain. The CR profile and spectrum depend on details of uncertain CR 
production mechanisms (including effects of the cluster assembly history), cooling processes, 
and transport processes via advection and diffusion [see, e.g., 74, and references therein]. 
In the case of nonrelativistic diffusive shock acceleration, the injection spectrum of CRs is 
expressed as s- m \ = (a + 2)/(<r — 1), where a is the compression ratio. A hard CR spectrum 
is achieved by high Mach number shocks that occur either at high redshifts during the 
formation of the protoclusters or today at the boundary where matter collapses from voids 
onto structured regions. AGN and supernova remnants may also inject CRs with rather hard 
spectra, which may be important [75], as well as the outflow from AGN [76]. Merger shocks 
with typical Mach numbers of M ~ 2 — 5 lead to steeper spectra. As noted above, the 
resulting CR spectrum inside GCs can be modified from the injection spectrum by various 
processes [36, 40]. 

In this work, for demonstration purposes, we simply use a power law, dS cr /dE p oc E~ s , 
where s is the CR spectral index. One should keep in mind that the CR spectrum can 
be made concave by detailed physical processes such as the nonlinear CR acceleration [e.g., 
33, 77]; but we avoid such model-dependent investigations. Then, the CR density n cr (r) is 
connected to the total CR energy as 



The spatial distribution of CRs is quite uncertain at present, so we consider two extreme 
limits with the CR spectrum independent of r. In one limit, expecting that CR transport 
is dominated by advection, we assume that CRs trace the thermal energy distribution. This 
would be the case especially for low-energy CRs. One numerical simulation showed that 
the cluster assembly history leads to a CR spectral index of s ~ 2.3 [40]. For simplicity, 
we consider the isobaric model, where the CR energy density is a fraction of the thermal 
energy density, X CT = U cr /Uth- Here U CT is the energy density of CRs above 1 GeV, whereas 
Uth = (3/2)[l+(l— 0.75xHe)/(l — 0.5xn e )]n e (r)kT e (r) is the thermal energy density [54, 55]. In 
another limit, active CR transport is relevant, and we assume that CRs are homogeneously 
distributed over the virial radius 2 . Such a situation may be realized at sufficiently high 
energies, especially in the > TeV range. For example, Ref. [36] suggested such a homogeneous 
distribution with s ~ 2.2 based on the radio and X-ray connection. In particular, in the 
scenario where CRs are accelerated up to VHE or UHE by accretion shocks [28, 43, 44], a 
CR population with a hard spectral index is expected, though weak shocks would dominate 
the gravitational energy dissipation in GCs. Note that the isobaric model can be regarded as 
an optimistic case in the sense that more gamma rays and neutrinos are produced compared 
to the uniform model. 

CRs are injected over a finite time, and spatial diffusion may lead to a break in the CR 
spectrum. When the cluster size is much larger than that of CR accelerators, in the steady 
state, diffusion can make the spectrum steeper by 1/3 for the Kolmogorov turbulence, i.e., 
s = Sinj + 1/3. In particular, if CRs diffuse without significant cooling after injections by a 

2 See Model B in Ref. [28]. Ref. [28] also considered the one-zone model by estimating the ICM density at 
the shock radius (Model A), but the external shock radius may be rather large, ~ 2r v i r [37, 42]. 




V 



(2.10) 
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central point source such as AGN in the CG core, one obtains n CT oc r" 1 if the steady state 
is achieved [26, 45]. In this case, depending on the CR energy, the upper limit of integration 
in Eq. (2.10) should be r max = min[rA c , r^is], where r^is is the diffusion length. We expect 
that the realistic situation would be between these two limiting cases. 

One purpose of this work is to demonstrate how we can constrain this parameter, i.e., the 
total CR energy, in a less model-dependent way through high-energy neutrino observations. 
The total CR energy can be related to the CR injection luminosity. For accretion/merger 
shocks, the accretion luminosity is L ac ~ 10 45 — 10 46 erg s _1 [35, 44], so the CR injection 
luminosity can be L cr ~ 10 45 erg s" 1 if ~ 10% goes to CRs. In addition, some AGN such as 
Fanaroff-Riley (FR) II galaxies, may have such high CR luminosities, and their luminosity 
could have been higher in the past. Supernova remnants may have a comparable energy 
budget, although CRs may suffer from adiabatic losses. In any case, the total CR energy 
stored in the CG is estimated to be 

£cr ~ L CI x min[i in j,t diff ] 

,61.5 / ^cv \ ( min[£ in j,£ diff ] 



where t; n j is the CR injection timescale, which is typically order of ~ 1 — 10 Gyr, and t<ns is 
the CR diffusion timescale. 

As in Eq. (2.3), the intensity (per unit energy, area, time and solid angle) is given by 

'* (e) = MiTTf J'Z d ' nMr) I iE '" (2 ' 12) 

where £ is the effective multiplicity of gamma rays or electrons/positrons or neutrinos, c is the 
relative velocity between an incident proton and a target nucleoli, and the CR normalization 

is given by f dE p (dS cv /dE p ) = 1. Here, Z max = yj r\ c — (dA@) 2 and = r/dA is the angular 
size measured from the cluster center. Examples of electromagnetic and neutrino spectra are 
shown in Figures 4 and 5. The calculation is performed based on SYBILL [78]. 

We evaluate the J factor for the CR case. As in the dark matter case, it depends on 
assumptions about the CR distribution. For CRs, the J factor defined by Eq. (2.8) has the 
following dependences, 



oc / dl n e (r), (2.13) 



27r0/" m 

9-ttP) T lso f 

oc J dl kT e (r)nl(r). (2.14) 

The results for the two CR profiles are shown in Figure 6. As shown, the best angular scale is 
model-dependent. For center-distributed CR profiles, the optimization should be done with 
smaller values of G. However, for the uniform CR distribution, the optimized angular scale 
is a few degrees. 

Though we focused on the pp reaction in this work, if the CR spectral index is hard 
enough and that the proton maximum energy is sufficiently high, the Bethe-Heitler process 
is important [44, 45]. Such limiting cases will be studied in future work. 

Also, one should keep in mind that gamma rays are produced via IC emission by elec- 
trons accelerated at accretion shocks [32, 34, 35]. Although the gamma-ray limit would give 
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Figure 4. Injection spectra of gamma rays (thick 
curves) and electron/positrons (thin curves) pro- 
duced by CRs via the pp reaction. The total CR 
energy is assumed to be £ cr = 10 61 erg for CR 
spectral indices of s = 2 and s = 2.25, and a uni- 
form CR distribution is used. 
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Figure 5. Injection spectra of neutrinos (sum 
of all flavors) produced by CRs via the pp re- 
action. The total CR energy is assumed to be 
£ cr = 10 61 erg for CR spectral indices of s = 2 
and s — 2.25, and a uniform CR distribution is 
used. 
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Figure 6. The J factor as a function of O, which represents the shape of the signal to the square root 
of the background at given energy, but now for the CR case. Two representative CR distributions are 
used, and the parameter set for Virgo is used. The typical angular resolution of IceCube (6 = 1 deg) 
and CTA (0 = 0.05 deg) are indicated by the dashed vertical line and the dotted vertical line, 
respectively. The dctectability is optimized at a few degrees for a uniform distribution. 



interesting constraints on the acceleration efficiency of electrons, it is beyond the scope of this 
work since our purpose is to demonstrate the power of neutrinos in constraining cosmic-ray 
ions in GCs. 
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2.3 Neutrino and gamma-ray emissions 

The formulas for calculating the injection spectra of neutrinos and gamma rays are given by 
Eqs. (2.3), (2.6) and (2.12). Before presenting results from dedicated numerical calculations, 
it is useful to give rough estimates. For annihilating dark matter, from Eqs. (2.4) and (2.6), 
one has (for gamma rays or neutrinos) 

F 2p _ l + b <OV >xx ( nnn n n M 20QC 2 



io n^r „™-2 / < av >xx \ m x c 



2.3 x 10" 10 GeV cm" 2 s" 



3 x 10~ 26 cm 3 s^ 1 J V 5 TeV 



,-2 



M 200 \{(l + b)g M \ fW\ /16Mpc x2 



5 x 10 14 M J V 10 4 J V n 



(2.15) 



where TZ = K(E) = (E 2 ^)' 1 f dE (Ej^) is the ratio of the bolometric flux to the differential 
flux at E (see Figures 1 and 2). 

For CRs, Eq. (2.12) leads to (for gamma rays or neutrinos) 



n N \f £ CT \ ( 20 \ /16 Mpc x 2 



~ 1.5 x 10" 1U GeV cm" 2 s" 1 ^ i: ^f i:r 12.16) 

V10" 4 cm- 3 / V10 61 er g y V^cr/ \ d / 

where 7?. cr = lZ cr (E p ) = (S 2 ^ 11 ) 1 J dE p {E p ^ L ) and we assume that the gamma ray or 
neutrino roughly carries a half the energy of the mesons. Since CRs are not depleted by the 
pp reaction, the pp meson production efficiency is estimated to be f pp (E p ) w nNKnPppCtmt, 
where K n a pp w 10 -26 cm 2 is the inclusive cross section of meson production, and t^t is the 
interaction time of CRs. If CRs are trapped in GCs, we expect t; n t ~ tmj, and the injection 
timescale may be the dynamical timescale of accretion/merger shocks, or the cluster age if 
CRs are continuously supplied by AGN inside the CG. If the CR escape becomes relevant, 
we expect t^t ~ rmnftdifi, r 20o/c]- In the latter case, when the diffusion coefficient increases 
with energy, gamma-ray and neutrino spectra become steeper than injection CR spectra. 

Dark matter and CRs may thus lead to comparable energy fluxes, even though both 
cases are affected by astrophysical uncertainties. Both predictions are not far from ranges of 
current and future experiments such as Fermi, CTA, and IceCube/KM3Net. 

Neutrinos reach Earth without attenuation. On the other hand, the fate of gamma 
rays is significantly affected by electromagnetic processes inside and outside the GCs. GCs 
typically have magnetic fields ~ 1 — 10 fiG around the cluster center, while ~ 0.1 — 1 /iG 
in the cluster outskirts [see, e.g., 79, 80, and references therein]. Since the Larmor radius of 
electrons/positrons is typically small enough, pairs produced in the GCs are trapped in the 
GCs, and lose their energies via the inverse Compton (IC) and synchrotron emissions 3 . Suf- 
ficiently high-energy gamma rays may also generate pairs via pair creation with low-energy 
photons, since the mean free path of gamma rays for the cosmic microwave background 
(CMB) is less than ~ Mpc at energies from ~ 100 TeV to ~ 3 EeV. As a result, IC emission 
is typically expected at E^ r w (4/3)7 2 £cmb — 3. 4 GeV (7 e /10 6 ) , where £cmb is the typical 



3 At UHE, the Larmor radius of pairs is rather large, and the synchrotron cooling length becomes much 
shorter. Then, the "synchrotron pair halo/echo" emission can be expected [81]. 
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CMB energy and 7 e is the electron/positron Lorentz factor. Then, attenuated and cascaded 
gamma rays leave the GCs. They may be further affected in intergalactic space. What hap- 
pens highly depends on uncertain intergalactic magnetic fields for whether the intergalactic 
cascade component contributes to the observed emission. To be conservative, we only take 
into account attenuation by the extragalactic background light (EBL) [see 82, and references 
therein] . 

In this work, the cascade gamma-ray spectrum inside the GC is numerically obtained 
by solving the following Boltzmann equations, 

-of = - i + ~it + -it + W ^ 

dN e _ dNT' on? a inj 

N * R K + -gf- - Q^[PsynN e ] + Q e ) (2-18) 



where 



f dn f du 
Rjj = / de— / — ca 77 (e,/i), 

Ric = J de^ J ^ cai C (e,/i), 



dt 
dNl c 
dt 



Here c = (1 — /x)c, t esc = T2oo/c is the photon escape rate, P syn is the synchrotron energy loss 
rate, iV 7 and N e are photon and electron/positron number densities per energy decade, and 
Q 1 ™ 1 and Q™ s are photon and electron/positron injection rates. For simplicity, we consider a 
one-zone model, assuming an average magnetic field of 0.3 (iG. The calculation is performed 
for the injection time of i; n j = 1 Gyr, and we obtain an essentially steady-state spectrum. 
For infrared and optical photons in GCs, we simply use the EBL with 10 times enhancement, 
and we adopt the low-IR model of Ref. [83]. Although this treatment is not fully accurate, 
it is sufficient for our purpose to demonstrate the power of neutrinos. 

Our results for the spectra from dark matter annihilation are shown in Figure 7 for 
cascaded gamma rays and Figure 8 for neutrinos. In the case of XX ~^ bb, the gamma-ray 
spectrum is dominated by gamma rays coming from hadronization, though there is a lower- 
energy tail due to IC cascades. For the XX ~ * channel, the IC emission gives the 
dominant contribution. In particular, CMB photons are upscattered by pairs coming from 
muon decay. 

Our results for the spectra from CRs are shown in Figures 9 for cascaded gamma rays 
and Figure 10 for neutrinos. The cascade component is relevant only when the CR spectrum 
is hard. For s = 2.25, cascade gamma rays are negligible in the GeV-TeV range, and they are 
relevant only below ~ 0.1 GeV. For s = 2, pionic gamma rays are still the most important, 
but the IC emission by pairs give moderate contributions, and its enhancement to the total 
gamma-ray flux is a factor of ~ 2 — 3. 
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Figure 7. Spectra of cascaded gamma rays 
produced by annihilating dark matter. The an- 
nihilation cross section is set to < av > xx = 
2.2 x 1CT 26 cm 3 s" 1 with substructures. The 
thick curves represent cases with EBL attenua- 
tion, while the thin curves are for cases without 
the EBL attenuation. For the bb channel, there 
is a high-energy bump due to pionic gamma rays 
and a lower-energy tail due to IC cascades. The 
Virgo cluster is considered, where the EBL at- 
tenuation is small. The Fermi limit is EFe ~ 
3 x 1CT 9 GeV cm' 2 s" 1 [51]. 



Figure 8. Spectra of neutrinos produced by 
annihilating dark matter. The annihilation cross 
section is set to < av > xx = 2.2 x 10~ 26 cm 3 s _1 
with substructures. The typical sensitivity of 
IceCube is EF E ~ 1(T 9 GeV cm" 2 s -1 . 




Figure 9. Spectra of cascaded gamma rays pro- 
duced by CRs via the pp reaction. We use the 
same injection spectra shown in Figure 4. The 
thick curves represent cases with EBL attenua- 
tion, while the thin curves are for cases with- 
out the EBL attenuation. The Fermi limit is 
EF E - 3 x 1CT 9 GeV enr 2 s" 1 [51]. 



Figure 10. Spectra of neutrinos produced by 
CRs via the pp reaction. We use the same injec- 
tion spectra shown in Figure 5. The typical sensi- 
tivity of IceCube is EF E ~ 1CT 9 GeV cm" 2 s" 1 . 



3 Neutrino and gamma-ray constraints 

In the previous section, we describe how we calculate neutrino and gamma-ray fluxes pro- 
duced by annihilating dark matter and CRs. Using these results, we derive neutrino con- 
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straints that could be placed by IceCube-like detectors in next several years, assuming the 
pessimistic case of non-detection. To obtain the neutrino constraints, we estimate the muon 
event rate by 

A^(> E v ) = SIT J dE' v A cS {E' u )^(E' u ), (3.1) 

where Q is the observed solid angle, T is the observation time, and A e s(E u ) is the neutrino 
effective area averaged over the zenith angle (that is not fully accurate in the UHE range). 
The full IceCube effective area is assumed to be three times as large as the IceCube-40 effective 
area [84, 85] with help of improvements in reconstruction techniques [58], and constraints 
are set by the criterion, iV S jg/ *J N s i g + iVbkg < 2. For the background iVbkg, we use the 
conventional atmospheric muon neutrino background [84, 86]. For the signal iV s jg, we consider 
neutrinos from annihilating dark matter or from CRs. Neutrino mixing is taken into account 
assuming 9±2 = 0.59 and 6*23 = vr/4 [87]. At low energies (at least in the < 300 TeV range), the 
atmospheric neutrino flux is dominant, and the constraints will improve as the square root of 
time. At higher energies, the atmospheric neutrino flux is negligible, and the constraints will 
improve linearly with time unless some other background is detected. The above criterion 
would be conservative at higher energies. 

We use T = 5 yr and = max[0 2 , 7r6^], where 6 U « 1.5y / TeV/£v deg is the kinematic 
angle. We take G = 0.5r2oo/^A; which is appropriate for annihilating dark matter with 
substructures and CRs with a uniform distribution (see Figures 3 and 6). Note that this 
gives conservative estimates for more center-distributed profiles. 

We only consider muon-track events. For cascade events, depending on the criteria 
of event cuts, the effective area is typically smaller than that for muon-track events [89]. 
Also, the angular resolution for cascade events is usually worse than that for muon-track 
events, which increases the atmospheric neutrino background though it can be lowered if 
muon neutrino charged-current events can be excluded [e.g., 90, 91]. In principle, the situ- 
ation can change if the angular resolution is improved in future neutrino detectors such as 
KM3Net (cf. Ref. [91]). Also, we do not consider any low-energy extension such as IceCube 
DeepCore [88], which can enhance the detectability below 100 GeV. 

Recently, gamma-ray constraints have been placed using Fermi and IACT data. Al- 
though we focus on neutrino constraints in this work, it is useful to compare. For that 
purpose, we show gamma-ray constraints taking into account cascades, from the flux upper 
limits used in Refs. [18, 51]. Although we do not perform dedicated analyses, this is enough 
for our present purpose of showing the neutrino power, and calculating the electromagnetic 
cascades allows us to see how the gamma-ray constraints extend to very high masses. 

In this work, we select five nearby clusters, Virgo, Fornax, Perseus, Coma and Ophi- 
uchus, using cluster parameters shown in Refs. [17, 54]. For cumulative neutrino backgrounds, 
see Refs. [28, 60]. We assume that no neutrino signal is observed in five years and consider 
either of the two contributions, annihilating dark matter or interacting CRs. In the future, 
some signal may be observed, and both may potentially contribute. However, the predictions 
of each contribution are currently uncertain and model-dependent, so that it is sufficient to 
focus on setting limits on either of them without any signal. If a signal is discovered, more 
detailed joint studies will be needed. 

3.1 Constraints on annihilating dark matter 

In Figure 11, we show forecasted neutrino constraints for five nearby clusters, assuming the 
channel. The Virgo cluster will give the most stringent bound, where the constraint 
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can be < av > X \ ^ 10~ 23,5 cm 3 s _1 at ~ 10 TeV in five years. For gamma-ray constraints, 
Virgo may not be a suitable target for the purpose of constraining dark matter because of 
contamination by gamma rays from the central AGN, M87. However, neutrino signals have 
not been observed yet, so we do not have to take care about the astrophysical background 
until it is detected. More distant clusters, Coma and Perseus, will give ~ 10 times weaker 
constraints, < av > xx < 10~ 22,5 cm 3 s -1 at ~ 10 TeV in five years. For Ophiuchus and 
Fornax in the southern hemisphere, we assume a hypothetical detector in the northern hemi- 
sphere. Our constraints agree with the independent work by Ref. [91] in the low-energy range, 
except that the neutrino effective area is slightly different and we consider the dependence of 
the atmospheric neutrino background on the zenith angle. There is a difference by a factor of 
~ 3 in the high-energy range, mainly because they took into account astrophysical neutrinos 
based on a model of Ref. [28] as a background. 

In Figure 12, we show neutrino constraints for five clusters for the bb channel. The 
constraints are much weaker than the results for the channel, because the spectrum 

due to hadronization is much softer than the spectrum due to lepton decay. Similarly to the 
channel, the Virgo cluster gives the most stringent limits, where the constraint can 
be < av >xx ~ 10 -22 cm 3 s _1 in five years. Note that, in Figures 15 and 16, the unitarity 
bound is also shown [92, 93]. Though this bound can be avoided if dark matter particles 
are not point-like, it typically gives the most stringent bound at sufficiently heavy masses, 
> 100 TeV. 

Our constraints can be compared to other neutrino constraints. The IceCube-59 anal- 
yses for dwarf galaxies gave < 10 -21 cm 3 s _1 for the t + t~ channel, which is weaker than 
the constraint for substructures but more model-independent [94]. The IceCube-22 analyses 
for the Galactic halo gave < 10~ 22 cm 3 s -1 for the [i~ channel and < 10 -20 cm 3 s" 1 for 
the bb channel, respectively [95]. The IceCube-40 analyses for the Galactic center obtained 
~ 10 times better constraints at ~ 1 TeV and comparable results at ~ 10 TeV [94]. In the 
presence of substructures with very small values of the minimum halo mass, the extragalactic 
component can be more relevant in diffuse backgrounds. For M mm = 1O _6 M0, the corre- 
sponding average flux-multiplier is go ~ (1 — 3) x 10 6 , and one has < 10 -23 cm 3 s" 1 for the 
channel and < 10 -21 ' 5 cm 3 s _1 for the bb channel, respectively, at ~ 30 TeV [60]. 

We here point out an interesting implication of our results. Recent preliminary Fermi 
data have indicated the possible existence of the "VHE Excess" in the ~ 100 GeV range [96] . 
This excess is attributed to cascaded gamma rays from astrophysical sources such as BL Lacs, 
but the dark matter interpretation could also be possible [60]. However, for annihilating dark 
matter, quite large cross sections are required. For go ~ 3 x 10 6 , one would roughly expect 
< av >x\ ~ 10 -22 ' 5 cm 3 s _1 for the p + fJ.~ channel and < av > xx ~ 10 -23 ' 5 cm 3 s _1 for the 
bb channel, respectively. However, the former case can be ruled out by neutrino observations 
of GCs, whereas the latter case is inconsistent with Fermi constraints from cascaded gamma 
rays (see below). However, note that accounting for the "VHE Excess" by dark matter decay 
is still allowed. 

In Figures 13 and 14, we show our results on individual clusters, Fornax and Virgo, 
respectively. The neutrino constraints are the same as those shown in Figures 11 and 12. 
We also show gamma-ray constraints (with cascades) estimated from Fermi data, although 
our limits are relatively conservative. Note that the latest gamma-ray constraints obtained 
by Refs. [19, 20] are better (though the constraints by Refs. [19, 20] is better than that by 
Ref. [21] partially due to different choices of A c ). Importantly, because the cascade effect is 
included in the calculation, our constraints extend to higher masses, even at > 10 TeV, and 



- 14 - 





-18 
-19 
-20 


'* "* """ ' ^ ^ 




-21 




CO 

E 


-22 


^^^^ — 


> 


-23 




O) 

o 


-24 
-25 


Virgo \. 

Fornax \. - 

Perseus \, 




-26 
-27 


Coma 

Ophiuchus 



2 3 4 5 6 

log(m [GeV]) 





- 1 o 

-19 


V I 






-20 






CO 

E 

U 


-21 
-22 






> 


-23 






_o 


-24 
-25 


Virgo 

Fornax 

Perseus 






-26 
-27 


Coma 

Ophiuchus '" , , , 



2 3 4 5 6 

log(m [GeV]) 



Figure 11. Forecasted neutrino constraints 
on the dark matter annihilation cross section, 
< av > xx , for five nearby GCs, with (thick 
curves) and without (thin curves) substructures. 
Nondctcction in five years is assumed for the 
channel. The unitarity bound (thick solid 
curve) and the typical cross section of thermal 
relics (thin solid curve) are also shown. 



Figure 12. The same as Figure 11, but for the bb 
channel. Only the constraints with substructures 
are shown since the constraints without them are 
quite weak. 



one sees that the constraints become weaker c< m x . 

The gamma-ray constraints are more stringent in the typically-considered range of the 
WIMP mass, < 100 TeV. In particular, for the bb channel, the neutrino bound is almost al- 
ways conservative. Even for the fJ- + ^~ channel, the gamma-ray constraints are more stringent 
below ~ 10 TeV, but the neutrino constraints are stronger at higher masses. This suggests 
that the neutrino constraints are the weakest at < 10 TeV, so neutrino observations typically 
give a conservative constraint on the total annihilation cross section [66, 92]. 

In Figure 14, for comparison, we also show the forecasted TeV gamma-ray constraint 
that can be placed by a future IACT, the Cherenkov Telescope Array (CTA) [97]. This 
indicates that the potential TeV gamma-ray constraints are also powerful, and it will provide 
better constraints for the bb channel compared to the neutrino constraint. Here, for simplicity, 
the sensitivity is reduced by 0/0 ps compared to the point source sensitivity, assuming the 
background-dominated case. We do not provide more detailed studies for such extended 
sources since the detailed performance of IACTs needs to be known, but it will be done in 
the future when one better knows the effective area, point spread function and backgrounds. 

Our forecasted neutrino constraints can be compared to other gamma-ray constraints. 
For example, Fermi constraints from dwarf galaxies at ~ 1 TeV are < 10 -22 ' 5 cm 3 s~ 1 for 
the channel and < 10 -24 cm 3 s _1 for the bb channel, respectively [98]. The con- 

straints become weaker cx; m x . The diffuse background measurements typically lead to 
weaker constraints [60, 99, 100]. Other limits from IACT observations also give interest- 
ing constraints [see a review 7, and references therein]. For the NFW profile, in the TeV 
range, the recent HESS measurement gave < 10 -21 ' 5 cm 3 s _1 for both the channels [49]. 
With significant boosts, < 10 -23 cm 3 s _1 for the bb channel was obtained. For the galactic 
center, a stronger constraint of < 3 x 10~ 25 cm 3 s" 1 is placed at TeV for the quark-antiquark 
channel [101]. Note that the isotropic DGB is difficult for IACTs to measure. 

Although we focused on the neutrino emission, the search for gamma-ray lines has been 
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Figure 13. Forecasted neutrino constraints 
(thick solid curves) and Fermi gamma-ray con- 
straints with cascades (thick dotted curves) on the 
dark matter annihilation cross section, < av > xx , 
for the Fornax cluster. The unitarity bound (thick 
solid curve) and the typical cross section of ther- 
mal relics (thin solid curve) are also shown. 
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Figure 14. The same as Figure 13, but for the 
Virgo cluster. In addition, we overlay future 
gamma-ray constraints that can be placed by 
CTA (thin dotted curves). 



interesting [102]. Recently, the existence of the line has been claimed around the Galactic 
center [103], and GCs may also be interesting sources [104]. 

3.2 Constraints on cosmic rays 

In Figure 15, we show the results of neutrino constraints on the total CR energy, f cr , assuming 
a uniform CR distribution. One sees that the neutrino bound is more stringent for harder 
CR spectral indices. This is simply because IceCube/KM3Net is most sensitive to ~ PeV 
neutrinos, where the atmospheric neutrino background is expected to be small. For spectral 
indices, s » 2, most of the CR energy is concentrated in the GeV range, so that the CR 
energy in the PeV range is so small that the neutrino constraints should be weak. One sees 
that the Virgo cluster gives the most stringent neutrino constraint, £ CT < 10 62 erg for s = 2. 

Neutrinos with ~ PeV energies are produced by protons with ~ 30 PeV [28]. Although 
it may be difficult to trap such high-energy CRs in GCs, it is useful to consider the isobaric 
model as an optimistic case. In this case, CRs are more clustered around the CG center, so 
the neutrino flux is enhanced for the same total CR energy. In Figure 16, we show forecasted 
neutrino constraints on the CR energy fraction in the isobaric model, X CT . More massive 
GCs are expected to be larger energy reservoirs and the neutrino flux is proportional to n 2 N 
rather than tin, so the order among the five clusters changes from that in Figure 15. One 
sees that the Perseus cluster gives the most stringent neutrino constraint, X CT < 0.03 for 
s = 2. 

In Figures 17 and 18, we show our results of multi-messenger constraints on CRs for 
Perseus and Virgo, respectively. The neutrino constraints are the same as those shown 
in Figures 15 and 16. To demonstrate the power of neutrinos, we also show gamma-ray 
constraints estimated from Fermi data, taking into account cascades. As expected, Fermi 
gamma-ray limits are sensitive to GeV gamma rays, so the constraints from gamma rays 
are more stringent for sufficiently steep CR spectral indices. For the Perseus cluster, the 
forecasted neutrino constraint is tighter than the gamma-ray constraint at s < 2.3. The 
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Figure 15. Forecasted neutrino constraints on 
the total CR energy, £ cr , for five nearby GCs. The 
uniform CR distribution is assumed. The Virgo 
cluster gives the most stringent constraint. The 
shaded region indicates the typical total CR re- 
quired in the scenario where GCs contribute to 
the observed CR flux. 
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Figure 16. Forecasted neutrino constraints on 
the CR energy fraction in the isobaric model, 
X CI , for five nearby GCs. The CR distribution is 
assumed to trace the thermal energy distribution. 
One sees that the Perseus cluster gives the most 
stringent constraint. 



relative importance seems weaker for the Virgo cluster, but we stress that the neutrino 
constraint provides an independent constraint on CR ions since GeV-TeV gamma rays can 
be produced by primary electrons and the extent of CR clustering may be different between 
GeV and PeV-EeV. Hence, independently of the gamma-ray limits, neutrino observations are 
important in order to constrain the amount of CRs that are accelerated by AGN or high- 
Mach number accretion shocks around the virial radius. In other words, our results suggest 
that the multi-messenger observations will be able to give interesting constraints on more 
specific models based on numerical simulations taking into account the cluster history, and 
details of CR acceleration and transport processes. 

In Figure 18, we show the forecasted TeV gamma-ray constraint that can be placed by 
CTA. It is indicated that the TeV gamma-ray constraint is also powerful, and it will provide 
better constraints for steeper CR spectral indices compared to the neutrino constraint. 

GCs could be sources of the observed CR flux. For CG accretion shocks, assum- 
ing a total mass accretion rate is M ~ O.lVg/G (where Vg is the velocity of the in- 
falling gas) [1], the dissipation rate through the accretion shock is estimated to be L ac ~ 

9x 10 45 erg s" 1 (/ 9 /0.16)M^ 3 , where f g = Qb/fl^ is the gas fraction. A fraction of the kinetic 
energy would dissipate by relatively high Mach-number shocks [37, 42]. For AGN, the typical 
jet luminosities of FR I and II galaxies are ~ 10 43 — 10 45 erg s _1 and ~ 10 45 — 10 47 erg s _1 , 
respectively [105]. In cither case, it is assumed that a significant fraction of the dissipated en- 
ergy is used for CR acceleration. On the other hand, the number density of GCs with masses 
above 1O 15 M is n cg sa 3 x 10~ 6 Mpc~ 3 [e.g., 106], but it becomes n cg sa a few xl0~ 5 Mpc -3 
for masses above 5 x 1O 14 M 4 . For GCs hosting AGN, only a fraction of GCs (and galaxy 
groups) would have powerful AGN, and n cg ~ 10 -5 Mpc -3 is used in Ref. [45]. Then, tak- 
ing into account the luminosity of CRs above 10 17 eV is smaller than that above GeV by 



4 Hence the prediction for individual GCs given by Ref. [28] is affected by the minimum mass of GCs that 
contribute to the observed CR flux, while the diffuse neutrino background prediction does not change much. 
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Figure 17. The forecasted neutrino constraint 
(thick solid curve) and the Fermi gamma-ray con- 
straint (with cascades) on the total CR energy, 
£ cr , for the Perseus cluster. One sees that the 
neutrino constraint is stronger for harder indices. 



Figure 18. The same as Figure 17, but for the 
Virgo cluster. For comparison, we overlay future 
gamma-ray constraints that can be placed by 
CTA (thin dotted curve). 



5 - 1000 (for s ~ 2 - 2.4), the energy budget of VHECRs may be 
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which can be comparable to the energy budget of observed CRs above ~ 10 eV, 
3 x 10 45 erg Mpc~ 3 yr -1 . Then, the diffuse neutrino flux can be order of > 
10~ 8 GeV cm" 2 , which could be seen by IceCube/KM3Net [28]. 

Next we briefly consider implications of future neutrino constraints on individual clus- 
ters. For example, for s = 2.25, the luminosity of injected CRs above above 10 17 eV, 
-^vhecr = 10 43 erg s" 1 , corresponds to L cr ~ 10 45 erg s _1 . For s = 2 below 10 17 eV and 
s = 2.5 above 10 17 eV [28], the corresponding luminosity becomes L CT « 10 erg s _1 . Then, 
through Eq. (2.11), the total CR energy amount may be £ CT fa 10 60 - 5 - 10 62 - 5 erg. Although 
details depend on the history of CR acceleration and escape properties, this implies that 
neutrino observations could test scenarios such that GCs contribute to the observed CR flux 
below the ankle. Note that only optimistic cases would be probed by IceCube/KM3Net via 
the search for individual steady sources, but stacking analyses can improve the situation. In 
addition, the diffuse flux limit would give useful constraints. 

In Figure 19, we show the case of lower and higher values of the proton maximum energy. 
For sufficiently low maximum energies, the constraint becomes weaker, since the atmospheric 
neutrino background is essentially critical. For higher maximum energies, the constraint does 
not change in the interesting range of the spectral index, s > 2, since the neutrino flux in the 
relevant energy range is almost the same. Note that, in this case, the constraint for s < 2 is 
optimistic due to severe attenuation in Earth. 

Finally, in Figure 20, we show the results for the CR energy fraction in the isobaric 
model, X CT , for the Virgo cluster. Again, independently of the gamma-ray limits, neutrino 
observations can provide constraints on the energy fraction of CRs that are accelerated by 
AGN or high-Mach number accretion shocks around the virial radius. 
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Figure 19. The dependence on the proton 
maximum energy in neutrino constraints for the 
Virgo cluster. The solid, dashed and dot-dashed 
curves are for 10 17 eV, 10 15 - 5 eV and 10 19 eV, 
respectively. 
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Figure 20. The forecasted neutrino constraint 
(thick solid curve) and the Fermi gamma-ray con- 
straint (with cascades) on the CR energy fraction 
in the isobaric model, X CI , for the Virgo clus- 
ter. We also overlay future gamma-ray constraints 
that can be placed by CTA (thin dotted curve). 



4 Summary and discussion 

In this work, we studied GCs as gigantic reservoirs of dark matter and CRs, and considered 
what we can learn from neutrino observations by IceCube and KM3Net, along with gamma- 
ray observations by Fermi and IACTs. Our work can be summarized as follows. 

(1) GCs could be interesting sources of neutrinos and gamma rays coming from an- 
nihilating dark matter, if the minimum halo size is very small. When the dark matter 
signal is boosted by substructures, the forecasted neutrino constraint for /U + /i~ channel gives 
< av > xx < 10 -23 cm 3 s _1 by five-year observations of IceCube or KM3Net, though it can be 
much weaker for other channels. We also studied the effects of electromagnetic cascades on 
the gamma-ray constraints, which allows us to compare between the neutrino and gamma-ray 
constraints at very high masses. For the bb channel, the neutrino constraint is almost always 
weaker than gamma-ray constraints. For channel, however, the neutrino constraint 
can be stronger than the gamma-ray constraints for m x c 2 > 10 TeV. At lower masses, where 
the gamma-ray constraints are stronger, the neutrino constraint should be regarded as a 
conservative bound on the total annihilation cross section. 

(2) GCs could be interesting sources of neutrinos and gamma rays since they are reser- 
voirs of CRs produced by accretion/merger shocks, hosted AGN and possibly supernova 
remnants. High-energy neutrino observations will give constraints on the total CR energy 
of £ cr < 10 62 erg with five-year observations of IceCube or KM3Net. If VHECRs trace the 
ICM gas, X cr < 0.1 can be obtained for s < 2.2, independently of gamma-ray observations. 
We also calculated the electromagnetic cascades, but they are typically irrelevant unless the 
CR spectrum is very hard. Although the Fermi gamma-ray constraint is stronger than the 
neutrino constraint for steep CR spectra, the neutrino limit is also powerful for CRs with 
hard spectra, which may be expected for CRs accelerated at accretion shocks or by AGN. 
Our study demonstrated that not only gamma-ray observations but also neutrino observa- 
tions give crucial constraints on the CR population in GCs, including the scenario where 
CRs from GCs explain the observed CR flux. Importantly, neutrino observations allow us 
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to probe not electrons but CR ions, unlike gamma-ray observations. Although we focus on 
the constraints assuming no signals are observed, possible detections in the near future may 
provide crucial clues to the nonthermal physics of GCs, the origin of CRs and so on. 

Multi-messenger observations, neutrino observations combined with gamma-ray, X-ray 
and radio measurements [107], will definitely enable us to give deeper insights into heavy 
dark matter and CRs contained in GCs. 

When we were completing this draft, a related but independent study by Dasgupta and 
Laha was posted on arXiv: 1206. 1322. 
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